Electronic response of aligned multishell carbon nanotubes 
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We report calculations of the effective electronic response of aligned multishell carbon nanotubes. 
A local graphite-like dielectric tensor is assigned to every point of the multishell tubules, and the 
effective transverse dielectric function of the composite is computed by solving Maxwell's equations. 
Calculations of both real and imaginary parts of the effective dielectric function are presented, for 
various values of the filling fraction and the ratio of the internal and external radii of hollow tubules. 
Our full calculations indicate that the experimentally measured macroscopic dielectric function of 
carbon nanotube materials is the result of a strong electromagnetic coupling between the tubes, 
which cannot be accounted for with the use of simplified effective medium theories. The presence of 
surface plasmons is investigated, and both optical absorption cross sections and energy-loss spectra 
of aligned tubules are calculated. 



The so-called multishell carbon nanotubes, with di- 
ameters on the order of 10 nm and lengths reaching up 
to several microns, are needle-like tubes consisting of a 
finite number of concentric cylindrical two-dimensional 
layers of graphite sheet arranged around a centraLhollow 
with a constant, separation between the layers.lxEl Since 
their discoveryjil these tubules have attracted much at- 
tention because of their sjpe. and their novel structural 
and electronic properties uu An important step in the 
investigation of these new materials was achieved when 
composites were manufactured with nanotubes organized 
into well-aligned arrays,Erlj since alignment of nanotubes 
was expected to be important not only to test proper- 
ties but also for potential applications such as atomic- 
scale field emitters. With the availability of aligned 
carbon nanotube films, optical—measurements were car- 
ried out with polarized light plil] thereby evaluating the 
frequency-dependent effective dielectric function of the 
composite and showing that carbon nanotubes have an 
intrinsic and anisotropic metallic behaviour. Accurate 
calculations of the effective dielectric function and the 
optical response of densely packed carbon nanotubes, as 
obtained by solving Maxwell's equations with the use 
of tensor- like dielectric functions, have been carried out 
only very recentlyJiil 

In this paper, we extend the calculations reported in 
Ref. [ll] to the case of hollow-core nanotubes, and em- 
phasize the effects of the inner radius on the effective 
dielectric function of the nanotube array. Calculations 
of the so-called energy-loss function, i.e., the imaginary 
part of the effective inverse dielectric function, are also 
presented, for various values of the filling fraction and 
the ratio of the internal and external radii of the tubules. 

We take a periodic array of infinitely long multishell 
nanotubes of inner and outer radii r and R, respec- 



tively, arranged in a square array with lattice constant 
a = 2 x R, as shown in Fig. 1. These tubules are assumed 
to be embedded in vacuum. In the energy range to 
10 eV of interest in the interpretation of absorption cross 
sections and energy-loss spectra, the diameter of typical 
multishell carbon nanotubes (2R ~ 10 nm) is small in 
comparison to the wavelength of light. Hence, our re- 
sults will be independent of the actual value of the outer 
radius of the cylinders. We also assume that this radius is 
large enough that a macroscopic dielectric function is as- 
cribable to the coaxial cylinders. These cylinders occupy 
a fraction / = 7r(l — j o 2 )/(4a; 2 ) of space, with p = r/R 
being the ratio between the inner and outer radii of the 
tubes. For simplicity, the magnetic permeabilities will be 
assumed to be equal to unity in all media. 

Planar graphite is a highly anisotropic material, and 
the dielectric function is a tensor. This tensor may be di- 
agonalized, by choosing cartesian coordinates with two of 
the axis lying in the basal plane and the third axis being 
the so-called c-axis. One defines the dielectric function 
s±(uj) perpendicular to the c-axis and the dielectric func- 
tion £|| (ui) for the electric field parallel to the c-axis. For 
carbon nanotubes, we assume full transferability of the 
dielectric tensor of planar graphite to the curved geome- 
try of carbon tubules, as suggested by Lucas et alLj for 
the case of multishell fullerenes. Hence, we simply assign 
a local graphite- like dielectric tensor to every point inside 
the nanotube and outside the inner core, 



£(oj) = e±(uj)(99 + zz) + £||(a>)rr, 



(1) 



99, zz, and rr being the unitary basis vectors of cylin- 
drical coordinates, and take the principal dielectric func- 
tions e±(u>) and £ii(u>) of graphite from Ref. [13[ 

In the long- wavelength limit, a composite material may 
be treated as if it were homogeneous, with the use of an 
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effective dielectric function e e //. The optical absorption 
cross section of the composite is then directly given by 
Im£ e //(w). Also, for small values of the adimensional 
parameter qR (qR < 1), q being the momentum transfer, 
the energy-loss spectra of a broad beam of swift electrons 
penetrating the composite is founcO to be well described 
by the q — » limit of the imaginary, part of the effec- 



tive longitudinal dielectric function,l£3 i.e., the so-called 
energy-loss function, Im[— e~ff(u)]. 

In order to compute, with full inclusion of the electro- 
magnetic interaction between the tubules, the effective 
dielectric function of our periodic system, we first solve 
Maxwell's equations. The local dielectric function of our 
composite material depends on frequency. T-hus, on-shell 
methods for solving Maxwell's equations,t3 which pro- 
ceed at constant frequency, are particularly well suited 
for our calculations. Though the methodology of Ref. [l6| 
was originally reported to solve Maxwell's equations with 
scalar dielectric functions, an extension of this formalism 
was developed in Ref. [ll], which enables us to consider 
tensor-like dielectric functions like the one of Eq. ([!]) . 

We first fix the frequency, thereby the local tensor-like 
dielectric function at any point of the composite material 
being specified, and approximate the continuous electro- 
magnetic field by its values at a series of mesh points 
located on a simple square lattice with unit cell of side 
a, as shown in Fig. 1. The corresponding discretized 
Maxwell's equations provide a relationship between the 
electromagnetic fields on either side of the unit cell, i.e., 
x and x + a (see Fig. I). We exploit this to calculate, 
on applying Bloch's theorem, the eigenvalues of the so- 
called transfer matrix, which will give the band struc- 
ture of the system, i.e., the dispersion relations k(ui) of 
the Bloch-state wave vector of propagation versus the 
frequency. The effective transverse dielectric function is 
then obtained as follows, 



k 2 {u)c 2 



(2) 



where c represents the speed of light. 

Fo an cm. wave normally incident on the structure 
(ky = k z = 0) there are two different values of e e ff(u) 
corresponding to s and p polarizations. For e.m. waves 
polarized along the tubes (s polarization), the electric 
field is parallel to the cylinders at every point, and is not 
modified by the presence of the interfaces. Hence, the 
effective dielectric function of the composite is simply 
the weighted-average of the dielectric functions of the 
constitucntsjia 



E 8 //(w) = /ex(w) + (l-/). 



(3) 



Our numerical results for the effective dielectric function, 
as obtained from Eq. (Q) for this polarization, accurately 
reproduce the exact results predicted by Eq. (J|), which 
represents a good check for our scheme. According to 
Eq. (H), Im£ e //(w) = /Imei(w); therefore, the shape 



of the experimentally determined Im£ e jj(w) must coin- 
cide for this polarization with Ira ex (a;). The measure- 
ments of e e ff(uj) reported in Ref. || for light polarized 
along the tubes are plotted in Fig. 2 by solid circles, to- 
gether with the effective dielectric function predicted by 
Eq. (||), which roughly reproduces with / ~ 0.5 the ex- 
perimentally determined lme e f f(uj). We note, however, 
that neither the exact position of the maximum nor the 
details of the real part of the measured effective dielectric 
function can be reproduced with the use of Eq. (||) . This 
may reflect existing differences between carbon structures 
in planar graphite and in nanotubes. Also, optical mea- 
surements may be affected by the presence of tubes with 
different orientations. 

The effective energy-loss function, lm[—e~h(u)], cor- 
responding to s polarization is also plotted in Fig. 2. For 
this polarization, the effective inverse dielectric function 
has a single pole, which occurs at a reduced plasma fre- 
quency where Ree_L ~ — (1 — /)// and Imei ~ 0. For 
graphite, at low concentrations of tubules (/ < 0.2) the 
resonance condition is never satisfied, and in the limit 
as / — > the effective energy-loss function simply coin- 
cides with ImEi(w), thereby showing the characteristic 
peak at ~ 4.6 eV. This peak is associated with the max- 
imum in the joint density of states of the n valence and 
conduction bands in graphite. 

Calculations of the p component of the effective di- 
electric function of an array of plain (p = 0) carbon 
nanotubes were reported in Ref. fl for various values 
of the ratio x between the lattice constant and the outer 
diameter of the cylinders, showing that the trend with 
increasing the concentration of tubules is for the actual 
dipolar peak in lme e ff(to) to be shifted from the isolated- 
cylinder dipole mode at ~ 6.5 eV to lower energies. When 
the nanoparticles are brought into close contact, elec- 
tromagnetic coupling between them converts the dipo- 
lar surface mode into a very localized one, trapped in 
the region between the nanostructures. At higher con- 
centrations of tubules, when graphite forms a connected 
medium (x < f), dipolar modes cannot be excited, and 
the optical absorption exhibits a single peak originated 
in the maximum of Ime(w) at ~ 4.6 eV. We have also 
calculated the enepgy-loss function of an array of plain 
carbon nanotubesj!3 and have found that for x < 1 it 
shows a single peak at the bulk plasmon resonance at 
~ 7eV, where Kee±(u>) and Im£i(w) are both small. 

For hollow tubes (p ^ 0), there are two distinct dipolar 
modes with either tangential or radial symmetry, simii 
lar to those present in the case of a thin planar filmcJ 
and a spherical shellEd With the presence of anisotropy 
this two-mode structure is replaced by a more compli- 
cated spectral representation^ Our full calculations of 
the p effective dielectric function of a periodic array of 
hollow carbon nanotubes, as obtained for various ratios 
between inner and outer radii, are shown in Figs. 3 
and 4 with x = 2.0 and x — 1.03, respectively^ As 
in the case of plain cylinders (p — 0), our full calcu- 
lations nearly coincide for x = 2.0 with a generalized 
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Maxwell-Garnett (MG) effective dielectric function ap- 
propriate for anisotropic hollow tubes.tfl However, in the 
close-packed regime (x = 1.03) the strong electromag- 
netic coupling between the tubes yields non-negligible 
contributions from multipolar resonances. This multipo- 
lar coupling provokes a redshift and a blueshift of the MG 
dipolar resonances that are visible in the optical spectra 
and the energy loss, respectively. These are the low- 
energy dipolar mode with tangential symmetry, which 
absorbs light (luie e ff(uj) is maximum at the correspond- 
ing energy), and the high-energy dipolar mode with ra- 
dial symmetry, which can be excited by moving charged 
particles (Im[-£~,,(w)] is maximum at the correspond- 
ing energy). 

The experimentally determined macroscopic dielectric 
function of close-packed carbon nanotubes, as reported 
in Ref. |9| for p-polarized light, is exhibited in Fig. 4 by 
solid circles. A comparison between these measurements 
and the p effective dielectric function of closed-packed 
(x = 1.03) plain (p — 0) carbon nanotubes was pre- 
sented in Ref. O, showing a nice agreement between the- 
ory and experiment, though a shift of approximately 1 
was needed to obtain agreement with the real part of 
e e //(o;)Ba If one chooses the ratio p between the inner 
and outer radii of the tubes to be 0.6, which yields in the 
close-packed regime (x — 1.03) a filling fraction / ~ 0.5, 
both real and imaginary parts of our calculated effective 
dielectric function show an excellent agreement with the 
experiment. This filling fraction (/ ~ 0.5) is precisely 
what is required to also reproduce the dielectric function 
experimentally determined with the use of light polarized 
along the cylinders, as discussed above (see Fig. 2). 

In conclusion, we have shown that the effective dielec- 
tric function of densely packed carbon nanotubes is the 
result of a strong electromagnetic coupling between the 
tubes. We have chosen the inner and outer diameters of 
hollow tubules to be 6 and 10 nm, respectively (p = 0.6), 
and have obtained both s and p effective dielectric func- 
tions to be in the close-packed regime (a = 10.3 nm, 
x = 1.03) in excellent agreement with the experimental 
measurements. Small discrepancies have been observed 
in the case of light polarized along the tubes (s polar- 
ization), which may reflect existing differences between 
the dielectric function of carbon structures in planar and 
curved graphite. 
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FIG. 2. The real 

and imaginary parts of the long-wavelength effective dielectric 
function, Ree e //(w) and Ime e //(w), and the energy-loss func- 
tion, 1xa.[—£~}f(u>)], of the periodic system described in Fig. 
1, for s polarized electromagnetic excitations. Dashed-dotted, 
long-dashed, short-dashed, dotted, and solid lines: calculated 
results, from either Eq. (||) or Eq. (0), for volume filling frac- 
tions of 19%, 35%, 47%, 65%, and 74%, respectively. The 
solid circles represent the experimental results reported in 
Ref. |. 

FIG. 3. The real 

and imaginary parts of the long-wavelength effective dielec- 
tric function, Ree e //(a>) and Ime e ;/(w), and the energy-loss 
function, Im[— e~* *(&)], of a periodic array of hollow carbon 
nanotubes, for p polarized electromagnetic excitations and for 
the ratio between the lattice constant and the outer diameter 
of the cylinders x — 2.0. Solid, long-dashed, short-dashed, 
and dashed-dotted lines represent our full calculations, as ob- 
tained from Eq. (0), for ratios between the inner and outer 
radii of the tubes p = 0.2, 0.4, 0.6, and 0.8, respectively. The 
dotted lines represent a generalized MG effective dielectric 
function appropriate for anisotropic hollow cylinders. 

FIG. 4. Same as Fig. 3, for the ratio between the lattice 
constant and the outer diameter of the cylinders x — 1.03. 
The solid circles represent the experimental results reported 
in Ref. §. 
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